Quantum criticality and the metal-insulator transition in 2D: a critical test. 
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Using recent insights obtained in heavy fermion physics on the thermodynamic singularity struc- 
ture associated with quantum phase transitions, we present here an experimental strategy to es- 
tablish if the zero-temperature transition in the disordered two dimensional gas is a real quantum 
phase transition. We derive a overcomplete set of scaling laws relating the density and temperature 
dependence of the chemical potential and the effective mass, which are in principle verifyable by 
experiment. 



The observation of a metal-insulator transition (MIT) 
in the two dimensional electron/hole gasses with 
quenched disorder (2DE/HG's) [1] revived the interest in 
the nature of disordered, interacting fermion systems [2] . 
A consensus has emerged that the transition is between a 
disordered Wigner crystal (Wigner glass) at low electron 
density and a strongly interacting Fermi-liquid state at 
higher densities. However, the nature of the phase tran- 
sition itself is subject of much debate. One viewpoint is 
that the interacting electron system near the transition 
is characterized by strong density inhomogeneities with 
puddles of high density electron gas immersed in a low 
density insulating background [3,4], suggesting a tran- 
sition of the percolation kind [5]. On the other hand, 
it has been suggested that the transition is a genuine 
continuous quantum phase transition (QPT) [6], finding 
support in for instance the scaling collapse of the resis- 
tivity [1,2]. In the past only transport quantities were 
accessible by experiment, but recently it appeared pos- 
sible to also measure thermodynamic quantities: the in- 
verse compressibility [7-10] and the effective mass [11,12], 
both of which appear to diverge near the transition. 

Thermodynamical properties are often easier to inter- 
pret than transport quantities. Very recently it was real- 
ized by Si, Rosch and coworkers in the context of heavy 
fermion physics that thermodynamics becomes extremely 
revealing in the vicinity of QPT's [13,14]. Based on al- 
most trivial scaling considerations, they demonstrated 
that an unambigous diagnostic can be constructed for 
a QPT, in the form of a plethora of scaling laws in- 
volving thermodynamic properties. The fact that ther- 
modynamics is more informative at a strongly interact- 
ing QPT than at a thermal phase transition is caused 
by the peculiar role of temperature as an effective fi- 
nite size in the former. The other important quantity is 
the zcro-tempcrature control parameter, tuning the sys- 
tem through the QPT. This coupling constant determines 
which thermodynamic quantities reveal the quantum sin- 
gularity. As we first realized in the context of cuprate 
superconductivity [15], the electron density is a partic- 
ularly interesting control parameter, because the singu- 



larity shows itself through the temperature- and density 
derivatives of the chemical potential, in addition to in- 
dependent information contained in the specific heat. It 
are precisely these quantities which have become accessi- 
ble by experiment in the MOSFET's and we will present 
here a recipy for how to use this information to uncover 
the QPT, if present in these systems. 

The thermodynamic singularity structure of a QPT is 
of a different kind than the one encountered at thermal 
phase transitions. The role of temperature is taken by 
the coupling constant y, a quantity tuning the system 
through its zero temperature transition at yc- We take 
as a working definition for a 'genuine' QPT that (a) at 
zero temperature the quantum dynamics \i\ D = d + z 
dimensional space time(2: is the dynamical exponent) be- 
comes scale invariant at y^ and (b) this critical state 
is universal in the sense that it obeys hyperscaling. In 
Euclidean path-integral representation, temperature sets 
the inverse conipactification radius of the imaginary time 
direction and acts therefore as a finite size. These pos- 
tulates can be compactly represented in the form of a 
scaling relation for the singular part of the free energy. 
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where r = (y — yc)/yc, the reduced coupling constant. 
Following Zhu et al. [13], Eq. (1) is equivalent to the 
following scaling forms for the free energy density. 
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where po and Tq are non-universal constants, while 
f{x) and f{x) are universal scaling functions. The first 
form is useful in the low temperature regime of the 
(dis)ordered phase in the proximity of the QPT, while the 
second form describes the quantum critical regime itself. 
Acoordingly, since there is no singularity at r = 0, T > 0, 
the 'quantum critical' function / can be expanded as 
fix ^ 0) ~ /(O) + xf{0) + (l/2)x2/"(0) + • • ■. On 
the other hand, in the stable phase the scaling function 
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is expanded as f{x) = /(O) + g{x) where g(x) describes 
the low temperature thermodynamics of the phases to 
the left- or right side of the QPT. By taking derivatives 
to the coupling constant and/or temperature one eas- 
ily derives the thermodynamic quantities. The fact that 
temperature enters through finite size scaling is a bless- 
ing; the outcome is a highly ovcrcomplctc set of scaling 
laws [13,15] relating the zero- and finite temperature ex- 
ponents of the thermodynamic quantities via yr and z. 

Let us now apply this general wisdom to the 2D 
MIT. The coupling constant clearly corresponds with 
the carrier density n while the transition occurs at a 
critical Uc such that the reduced coupling constant is 
r = (n — nc)/nc. The thermodynamic quantities as- 
sociated with the density are the chemical potential 
/i = dF / dn and the inverse compressibility (or incom- 
pressibility) (I/k) = {l/n){d'^ F / dn?). In addition one 
has the specific heat C = -T{d'^F/dT'^) and the tem- 
perature derivative of /i, T„ = (dfi/dT) = (d'^F/dTdn). 
In principle, one can derive a total of 9 scaling forms 
for the singular parts of 1/k, t] and C, applicable to the 
quantum critical regime (r = 0, finite T) and the low 
temperature regimes (T ^ 0, r 7^ 0) of the two phases 
in the proximity of the QPT, respectively. Although the 
physics is in any other regard very different, the 2d MIT 
shares the density as coupling constant with the (sus- 
pected) QPT in optimally doped cuprates [16] and one 
can find a full set of scaling forms in a paper [15] dealing 
with the specifics of the cuprates. 

Another difference with the cuprates is that the quanti- 
ties of crucial importance (compressibility, chemical po- 
tential) are now routinely measured in the MOSFET's 
[7-10] while this information is not available yet in the 
cuprates [15]. In the remainder we will be primairely 
interested in the high density side of the MIT. All we 
have to know is the function g{x) in this regime, and 
this can be deduced from the low temperature behav- 
ior of the specific heat C. Because of their 'lack of 
volume' C cannot be measured directly and we will as- 
sume that this state is 'like a Fermi-liquid' in the sense 
that the following Fermi-liquid relation is still valid: 
(C/Co) = 1 + Ff/S = (mVm), where Co = vT, the 
specific heat in the non-interacting limit is the density 
of states) while Ff is the Landau parameter which also 
governs the transport effective mass m* /m. This effec- 
tive mass is directly measured in the magneto-transport 
experiments [12], and found to bo finite at zero temper- 
ature, suggesting a linear specific heat C — jT. It is 
easily checked that this implies g{x) = ca:^°+^ (masless) 
spectrum with uq = 1. 

The scaling relations for the various quantities are eas- 
ily obtained and let us directly specialize to the case of 
the 'Fermi-liquid like' regime in the proximity of the 
2D MIT {d = 2 and yo = 1). The unknown expo- 
nents are the dynamical exponent z and the coupling 
constant (electron density) exponent yr- It is convenient 



to reparametrize the latter in terms of an exponent hav- 
ing a similar status as the specific heat exponent a in the 
case of a thermal phase transition [15], 
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The singular, low temperature contributions to the 
various thermodynamics quantities in the liquid have the 
following form. 
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where the various coefficients which depend critically 
on r = {n — nc)/nc are given by, 
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One infers that the divergence of the prefactors 
7, u, Kt of the temperature dependent parts of C, T and 
1/k are all governed by the fundamental mass exponent 
defined through ^ m* /m = r"'" with 



am = {2-ar){2-z)/{2 + z), 



such that V 



and Kt 



(6) 

These rela- 



tions between the exponents are just a consequence of 
the assumption that r enters the free energy in the form 
of a power law and the Fermi-liquid form for /. Hence, 
by measuring with high precision the temperature depen- 
dence of the chemical potential and incompressibility one 
can critically test if the system is like a Fermi-liquid with 
a mass diverging at the MIT. According to the transport 
measurements [11,12] m* ~ r"" with a™ — —0.5 ±0.1 
and when the interpretations are correct it has to be that 

In fact, the first genuine scaling laws in this liquid 
regime can be deduced from universal amplitude ratio 's. 
Usually amplitudes are non-universal but Zhu et al. [13] 
realized that in the case of pressure driven QPT's the 
Grueneisen ratio becomes universal. The analogon of the 
Grueneisen ratio in the case of the density driven transi- 
tion is the ratio of T and C which is expressed entirely in 
the exponents while the non- universal factors (^05^0,0) 
drop out. 



{2^z){2-ar)\ _i 
1 J. 



(7) 



2 



In addition, Eq. (5) implies a second universal ratio 
which has not been identified before, 



Tt; ^ / 2(2 -a.) 

Kt \{2-ar){2-z)-2-z 
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It follows from the scaling analysis that the divergence 
of the zoro-tcmpcraturc incompressibility Kq is governed 
directly by the elementary 'energy exponent' - this 
is no wonder since the incompressibility has the same 
status as a specific heat when the roles of temperature 
and density are exchanged. Hence, the zero temperature 
incompressibility is the quantity which reveals the quan- 
tum singularity most directly in case of density-driven 
transitions. As an interesting aside, the screening length 
s ~ I/k in a 2D charged liquid and this length is therefore 
governed at zero temperature by the exponent a^. One 
could be tempted to identify the screening length with 
the correlation length but this is not quite right: the 
correlation length is governed instead by the exponent 
v = l/Ur- In fact, the correlation length is the length 
scale associated with order parameter correlations and 
the screening length has a completely different, 'truely 
thermodynamic' status. There might well be an order 
parameter at work at the MIT but we have left it com- 
pletely implicit in the present analysis which is equivalent 
to the finite size scaling of the specific heat singularity at 
a thermal phase transition. 

Up to this point we have indentified three scaling laws 
governing the liquid: one relating a^. to the mass ex- 
ponent am via Eq. (6), and the universal ratio's Eq.'s 
(7,8). This is not all, because an other set of scaling re- 
lations can be deduced governing the quantum critical 
regime itself. At the quantum-critical point the scales 
associated with the stable fixed points should collapse. 
Approaching the MIT from the metallic side the situa- 
tion is in this regard analogous to what is found in the 
heavy fermion systems: the Fermi-energy (the scale of the 
Fermi-liquid) has to vanish [17], or equivalently, the effec- 
tive mass has to diverge. In the quantum critical regime 
the temperature dependence of the thermodynamic prop- 
erties are governed by a different set of scaling laws which 
are also governed by the fundamental exponents {z,ar)- 
The quantum critical regime can be accessed right at the 
QPT, but also away from the critical point, as long as 
temperature exceeds the Fermi-energy [6]. Hence, if the 
MIT is governed by a genuine QPT it has to be that upon 
approaching the phase transition a temperature window 
opens up of universal behaviors bounded from the below 
by the Fermi-energy and from the above by a cross-over 
to a non-universal ultraviolet. 

The temperature dependences of C, T and I/k in this 
quantum critical regime are easily deduced from the sec- 
ond form for the singular part of the free-energy in Eq. 
(2), 
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It follows that the temperature exponent of the spe- 
cific heat directly reveals the dynamical exponent z - a 
direct consequence of the role of temperature as a finite 
size. Since there is no a-priori relation between trans- 
port and specific heat in the quantum critical regime it 
appears that this quantity is beyond the reach of exper- 
iment. However, again the chemical potential is a valu- 
able source of information. Its temperature and density 
derivatives should obey two other independent scaling 
laws. 

In summary, we have identified five scaling relations 

depending on two unknown exponents {ar,z), govern- 
ing the behaviors of quantities which can be measured in 
principle by present day experimental techniques, which 
should be obeyed if a genuine QPT is at work in the 
2DEG's. How does the above relate to the thermody- 
namic experiments, feasible in the MOSFET's? High pre- 
cision measurements of the effective mass and the chem- 
ical potential in the same samples are required, to es- 
tablish if these quantities reveal divergences at the same 
critical density. The dynamical range should be improved 
to densities much closer to the transition and a main 
message of the present work is that the temperature de- 
pendence of especially the chemical potential contains 
valuable information. 

What can be said on basis of the data which are at 
present available? The divergence of the effective mass 
seems to be reasonably well established, suggesting a 
mass exponent am = —0.5 ± 0.1 [11,12]. The incom- 
pressibility data suggest that Kq diverges at some crit- 
ical density [7 10]. However, the dynamical range cov- 
ering this divergence is in the present data too small. 
Pending the choice of Uc we find that the incompress- 
ibility upturn can equally well be fitted with an in 
the range 0.2 — 2 [18]. As we already emphasized, is 
analogous to the specific heat exponent of thermal phase 
transitions, and these are generically small compared to 
unity. Even within the uncertainties, the present data 
seem to suggest that the incompressibility divergence is 
much stronger than the specific heat divergence at a typ- 
ical classical phase transition. It is easy to derive some 
bounds on the allowed values of the exponents - if these 
are violated one is surely not dealing with a QPT. The 
fundamental exponent is the coupling constant dimen- 
sion yr = {d+ z)/{2 — ar). To ensure the relevancy of 
the coupling constant yr > and this implies an upper 
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bound a,. < 2. 

Another bound follows from the condition that both 
the effective mass and the incompressibility are diverg- 
ing at the transition. The divergence of the mass implies 
am = —\ctm,\j while a diverging incompressibility needs 
Ur > 0. Using the scaling law Eq. (6) to express ar 
in terms of a^n as a,. = 2 + jam] (2 + z)/(2 — z), we 
find a lower bound on z > (2 + |am|)/(l — |ar„|/2). For 
instance, taking the experimental value \am\ — 1/2, it 
follows that z > 10/3. This is a surprisingly large lower 
bound, once again demonstrating that quantum phase 
transitions involving fermions have little in common with 
their bosonic- or classical sibblings [17]. 

Summarizing, we have presented in the above a pow- 
erful diagnostic to establish whether the metal-insulator 
transition in 2DEG's is a genuine quantum phase tran- 
sition. We are not advocating the view that it has to be 
a quantum phase transition. The only way to find it out 
is by more and better experiments. Alternatively, one 
could take the viewpoint that the MIT cannot possibly 
be a QPT given the intrinsic tendencies of the system 
to micro-phase separate. The various arguments [3,4] 
in support of this viewpoint are persuasive, but only so 
to the extent that the starting assumption is satisfied: a 
good metal separated from an insulating crystal by a first 
order transition which is in turn 'masked' by electrostatic 
effects and/or the disorder potential. In this context, the 
'quality' of the metal refers to its screening length being 
small. In 2D, the screening length s ~ 1/k and accord- 
ing to a Fermi-liquid relation v/k = (m/m*)(l -f Fq). 
Even when Fq can be ignored, the incompressibility and 
thereby the screening length has to diverge in the ap- 
proach to the transition, because the effective mass is di- 
vergent. It is hard to imagine how a mass divergence can 
arise from inhomogeneities - it is most likely associated 
with the physics of the uniform liquid [19]. Stronger, 
as is obvious from the scaling analysis, the screening 
length and the mass are in principle governed by differ- 
ent exponents, which implies that also Fq might behave 
singularly. Microscopically this might have to do with 
the effect discovered by Das Sarma in the 1980's that 
quenched disorder enhances the screening length [20], a 
physics which was recently discussed in the context of 
the MIT by Si and Varma [21]. The bottomline is that 
the homogeneous liquid has an intrinsic tendency to be- 
come increasingly incompressible approaching the MIT, 
and thereby increasingly similar to the Wigner glass in its 
electrostatic responses. This in turn will tend to diminish 
the tendency towards the formation of inhomogeneities, 
potentially paving the road for a continuous quantum 
phase transition. In an upcoming publication we will ad- 



dress these matters on a more quantitative level. 
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